Transgenerational Paternal Inheritance of TaCKX GFMs Expression Patterns Indicate a Way to Select Wheat Lines with Better Parameters for Yield-Related Traits

Members of the TaCKX gene family (GFMs) encode the cytokinin oxygenase/dehydrogenase enzyme (CKX), which irreversibly degrades cytokinins in the organs of wheat plants; therefore, these genes perform a key role in the regulation of yield-related traits. The purpose of the investigation was to determine how expression patterns of these genes, together with the transcription factor-encoding gene TaNAC2-5A, and yield-related traits are inherited to apply this knowledge to speed up breeding processes. The traits were tested in 7 days after pollination (DAP) spikes and seedling roots of maternal and paternal parents and their F2 progeny. The expression levels of most of them and the yield were inherited in F2 from the paternal parent. Some pairs or groups of genes cooperated, and some showed opposite functions. Models of up- or down-regulation of TaCKX GFMs and TaNAC2-5A in low-yielding maternal plants crossed with higher-yielding paternal plants and their high-yielding F2 progeny reproduced gene expression and yield of the paternal parent. The correlation coefficients between TaCKX GFMs, TaNAC2-5A, and yield-related traits in high-yielding F2 progeny indicated which of these genes were specifically correlated with individual yield-related traits. The most common was expressed in 7 DAP spikes TaCKX2.1, which positively correlated with grain number, grain yield, spike number, and spike length, and seedling root mass. The expression levels of TaCKX1 or TaNAC2-5A in the seedling roots were negatively correlated with these traits. In contrast, the thousand grain weight (TGW) was negatively regulated by TaCKX2.2.2, TaCKX2.1, and TaCKX10 in 7 DAP spikes but positively correlated with TaCKX10 and TaNAC2-5A in seedling roots. Transmission of TaCKX GFMs and TaNAC2-5A expression patterns and yield-related traits from parents to the F2 generation indicate their paternal imprinting. These newly shown data of nonmendelian epigenetic inheritance shed new light on crossing strategies to obtain a high-yielding F2 generation.


Introduction
Bread wheat (Triticum aestivum) is the most important cereal crop in the temperate climate and provides a staple food for more than a third of the world's population [1]. It belongs to the Triticeae tribe, which also includes barley, rye, and triticale. Among the species, wheat has the largest and most complex hexaploid genome (2n = 6x = 42), which consists of the three homoeologous subgenomes A, B, and D. Each gene present as homologues A, B, and D could retain its original function or, as a result of independent evolution, develop heterogeneous expression, and/or one or two copies may be silenced or deleted [2,3].

Cooperating and Opposite-Functioning Genes
TaCKX5 with TaCKX9 (yellow) and TaCKX2.1 with TaCKX11 (green) showed coordinated up-or downregulation in 7 DAP spikes of the paternal parent of C1, C2, C6, and C8 crosses; and C3, C4, C7, and C8 crosses, respectively ( Table 2). Higher expression of TaCKX5 and 9 in this parent was associated with a higher yield in F2. However, a higher coordinated expression of TaCKX2.1 with TaCKX11 in the paternal parent determined a lower yield in F2, and, in contrast, a lower expression of these two genes in the paternal parent was associated with a higher yield.
In the 7 DAP spikes, the paternal parent of the C3, C4, C5, and C6 crosses, TaCKX2.2.2, showed opposite expression to TaCKX10 (blue), and upregulation of the first and downregulation of the second were associated with lower yield (but not in C6). In the paternal parent of the C1, C2, C7, and C8 crosses, NAC2 was oppositely expressed to TaCKX9; in these crosses, a high yield was observed when TaCKX9 was upregulated and NAC2 was downregulated, and vice versa (only in C7 and C8). Furthermore, upregulated TaCKX5 and downregulated TaCKX1 were associated with high root mass in C2 and conversely in the reverse cross (C1).
Among the TaCKX genes coordinately expressed in the paternal seedling roots were TaCKX3, 5, and 8 (green) in C5, C6, C7, and C8 crosses; TaCKX3 and 8 (green) in C1 and C2 crosses; TaCKX10, 11, and 1 (yellow), and NAC2 in C3 to C7 crosses; and TaCKX10 and 11 (yellow) in C1 and C2 crosses. However, in one reciprocal cross, C3 and C4, the expression of TaCKX3 and TaCKX5 was opposite, and in the case of upregulation of TaCKX3 and downregulation of TaCKX5, the grain yield in F2 was higher. Depending on the crosses, downregulation of TaCKX5 with TaCKX9 and upregulation of NAC2 in spikes of the low-yielding maternal parent and the opposite regulation of these genes in spikes of the higher-yielding paternal parent resulted in high-yielding F 2 , characterized, as in the paternal component, by a higher expression level of TaCKX9 and a lower expression level of NAC2. Upregulation of TaCKX2.1 and 11 in spikes of the maternal parent and downregulation of these genes in the paternal parent were associated with downregulation of TaCKX11 in high-yielding F 2 . Similarly, the upregulation of TaCKX2.2.2 and the downregulation of TaCKX10 in the spikes of the low-yielding maternal parent and opposite regulation of these genes, and the yield in the paternal parent, resulted in the downregulation of TaCKX2.2.2 and the upregulation of TaCKX10 in the spikes of the high-yielding F 2 . The expression of TaCKX3 in seedling roots of the high-yielding paternal parent and F 2 was upregulated. However, TaCKX8 expression was upregulated, and NAC2 was downregulated in the same organ of the paternal parent, but these genes were up-or down-regulated in F 2 , depending on the cross.

Cooperating and Opposite-Functioning Genes
TaCKX5 with TaCKX9 (yellow) and TaCKX2.1 with TaCKX11 (green) showed coordinated up-or downregulation in 7 DAP spikes of the paternal parent of C1, C2, C6, and C8 crosses; and C3, C4, C7, and C8 crosses, respectively ( Table 2). Higher expression of TaCKX5 and 9 in this parent was associated with a higher yield in F 2 . However, a higher coordinated expression of TaCKX2.1 with TaCKX11 in the paternal parent determined a lower yield in F 2 , and, in contrast, a lower expression of these two genes in the paternal parent was associated with a higher yield.
In the 7 DAP spikes, the paternal parent of the C3, C4, C5, and C6 crosses, TaCKX2.2.2, showed opposite expression to TaCKX10 (blue), and upregulation of the first and downregulation of the second were associated with lower yield (but not in C6). In the paternal parent of the C1, C2, C7, and C8 crosses, NAC2 was oppositely expressed to TaCKX9; in these crosses, a high yield was observed when TaCKX9 was upregulated and NAC2 was downregulated, and vice versa (only in C7 and C8). Furthermore, upregulated TaCKX5 and downregulated TaCKX1 were associated with high root mass in C2 and conversely in the reverse cross (C1). Among the TaCKX genes coordinately expressed in the paternal seedling roots were TaCKX3, 5, and 8 (green) in C5, C6, C7, and C8 crosses; TaCKX3 and 8 (green) in C1 and C2 crosses; TaCKX10, 11, and 1 (yellow), and NAC2 in C3 to C7 crosses; and TaCKX10 and 11 (yellow) in C1 and C2 crosses. However, in one reciprocal cross, C3 and C4, the expression of TaCKX3 and TaCKX5 was opposite, and in the case of upregulation of TaCKX3 and downregulation of TaCKX5, the grain yield in F 2 was higher.
In three reciprocal crosses, NAC2 was downregulated in paternal roots (C2, C4, and C8), and in two of them (C2 and C8), NAC2 was downregulated in paternal spikes as well. This negative regulation of NAC2 occurred in the F 2 progeny, which was accompanied by a higher yield and a higher or similar to the parents' mass of the seedling roots. In contrast, in another way crosses, when expression of NAC2 was increased in paternal roots (C1, C3, and C7) and was upregulated in paternal spikes, the same was observed in F 2 progeny, characterized by lower yield and lower or similar to the parent mass of the roots.
The higher yield in F 2 has been associated with the same or higher CKX activity, as in the paternal parent, in 7 DAP spikes. A higher number of semi-empty spikes, which occurred in low-yielding F 2 of the C3, C5, and C7 crosses, was accompanied by downregulated TaCKX10 and/or upregulated TaCKX11 in 7 DAP spikes and upregulated TaCKX10, 11, and NAC2 or downregulated TaCKX10 and NAC2 in seedling roots.

The Correlation Coefficients between TaCKX GFMs and NAC2 Expression, CKX Activity, and Yield-Related Traits Were Significant for Both Parents or the Maternal or Paternal Parent Separately
The correlation coefficients between TaCKX GFMs and NAC2 expression, CKX activity, and yield-related traits in reciprocal crosses were analyzed separately for the maternal parent and F 2 , paternal parent, and F 2 for each cross (Table S1).

Correlations between TaCKX GFM and NAC2 Expression and Yield-Related Traits in the Group of Maternal Plants, and F 2 and Paternal Plants, and F 2 of Reciprocal Crosses
Seed number and spike number were positively correlated (Tables 3 and S1); however, each of these yield-related traits was correlated with different TaCKX GFMs.    All correlation coefficients ≥0.60; bold-significant correlation coefficients; nc-no correlation; +-positive correlation; =+-low positive correlation; ++-strong positive correlation; −-negative correlation; =−-low negative correlation; −−-strong negative correlation.

Seed Number
The decrease in seed number in maternal plants and their F 2 (M and F 2 ) of the C1 cross (S12B × S6C) was strongly negatively correlated with upregulated TaCKX1 and TaCKX5 in spikes and positively correlated with the downregulated TaCKX11 in the seedling roots of F 2 . There was no significant correlation between the expression of TaCKX GFM and the yield-related traits in the groups of paternal plants (P) and F 2 in the same cross, and M and F 2 , and P and F 2 in the reverse, C2 cross. The F 2 progeny in this reverse cross showed a similar yield and greater root mass compared to the parents.
In both M and F 2 , and P and F 2 of C3, the decrease in seed number was strongly positively correlated with TaCKX2.1 and TaCKX2.2.2 in spikes and positively correlated with TaCKX3 and TaCKX8 but negatively correlated with TaCKX11 in seedling roots. These correlations were not significant in the reverse, C4 cross, in which the M (KOH7) and F 2 plants showed higher yields.
The decrease in seed number in M × F 2 of C5 was negatively correlated with TaCKX1 in spikes and positively correlated with downregulated TaCKX3 and 5, and upregulated NAC2 in seedling roots. There were also positive correlations of TaCKX3 in roots between P and F 2 of the same cross. These correlations were not significant in the reverse C6 cross; however, F 2 of this cross was characterized by higher yield and similar root mass than in the parents.
The increase in seed number was strongly positively correlated with downregulated TaCKX2.1 in spikes and negatively correlated with downregulated TaCKX1 in the seedling roots only in F 2 progeny of a C8 cross, 14K (in the case of P × F 2 only for TaCKX2.1).

Spike Number
The decrease in the number of spikes in M × F 2 and P × F 2 of the C1 cross was strongly positively correlated with TaCKX2.1 and TaCKX2.2.2 in spikes and positively correlated with TaCKX11 in seedling roots (only in M × F 2 ). There was also a significant and positive correlation between the expression of TaCKX2.2.2 and the number of spikes in the reverse C2 cross of M and F 2 . Furthermore, in the same cross, the spike number was negatively correlated with downregulated NAC2, but only in the P × F 2 group. There were no significant correlations between TaCKX GFM expression and spike number in M × F 2 and P × F 2 spikes of C3 and C4. However, there was a strong and positive correlation of the spike number with TaCKX8 in the roots of C3 and TaCKX5 in the roots of C4.
The decreased spike number in M × F 2 and P × F 2 of C5 was not correlated with any TaCKX expressed in the spikes but was negatively correlated with TaCKX1, positively correlated with TaCKX5, and positively correlated with NAC2 expressed in the seedling roots. Conversely, in reverse C6 cross, there was a positive correlation of the spike number with TaCKX2.1 in a P × F 2 , which resulted in a higher yield phenotype in the F 2 .
The spike number in C7 and C8 crosses was not correlated with the level of expression of any gene tested in the spikes; however, it was negatively correlated with the expression of TaCKX1, 5, and 8 in seedling roots.

TGW
TGW was positively correlated with TaCKX2.1, 10, and NAC2 in spikes of M and F 2 of C1 and with TaCKX2.2.2 of P and F 2 of the same cross. There was no correlation in F 2 between the expression of the genes tested and TGW in spikes of the C2 and roots of the C1 and C2 crosses. There were no correlations between the TGW and TaCKX genes in the spikes and roots of C3. However, there was a negative correlation of this trait with TaCKX2.2.2 in spikes and positive correlations with TaCKX10 and NAC2 in roots of the reciprocal C4 cross. Positive correlations of NAC2 with TGW were also observed in the roots of C5 but not in those of C6. Negative correlations of TaCKX2.1 and 2.2.2 in spikes with TGW were observed in both reciprocal crosses, C7 and C8. Additionally, TaCKX9 was negatively correlated with the trait in M + F 2, TaCKX5 was positively correlated with the trait in P + F 2 of C7, and TaCKX10 was negatively correlated with TGW in P + F 2 of C8. There was no correlation between TGW and any gene expression in roots.

Root Mass
The mass is positively correlated with the expression of TaCKX5, 11, and NAC2 in spikes of M + F 2 of C1 and negatively correlated with NAC2 in M + F 2 of C2. A positive correlation between root mass and TaCKX11, and NAC2 was also visible in P + F 2 of C3, and a negative correlation between trait and NAC2 expression was also observed in spikes of M + F 2 of C5 and P + F 2 of C7. Furthermore, in the C1 cross, this trait was negatively correlated with the expression of TaCKX1, 3, 10, and NAC2 in roots of M + F 2 and with TaCKX11 in roots of P + F 2 . There were also positive correlations between root mass and TaCKX1 (M + P of C3), root mass and TaCKX2.1 (P + F 2 of C4; P + F 2 of C 5 ; M + F 2 of C6), and root mass and TaCKX2.2.2 (P + F 2 of C4; M + F 2 and P + F 2 of C6). The expression of another gene, TaCKX10, in spikes, was positively correlated with root mass in P + F 2 of C4 but negatively correlated in M + F 2 of C8. Correlations between root mass and gene expression tested in roots were dependent on the parent and cross. There were negative correlations with TaCKX1 in 5 out of 16 combinations tested, negative correlations with TaCKX3 in 3 combinations, but positive correlations in two combinations, positive correlations with TaCKX5 in two combinations, and single positive or negative correlations with TaCKX8, 11, and NAC2. The root mass in single combinations was positively correlated with the yield (twice), height of the plant (once), and length of the spike (once), and negatively correlated with seed number of seeds (twice).

Semi-Empty Spikes
The number of semi-empty spikes was positively correlated with the expression of TaCKX9 in the P + F 2 C1, C2, and M + F 2 C3 crosses, positively correlated with TaCKX5 in the M + F 2 , and P + F 2 C1 and C6 crosses, and positively correlated with TaCKX10 in the M + F 2 C7 and P + F 2 C8 crosses, all expressed in 7 DAP spikes. The negative correlation between the number of semi-empty spikes and TaCKX2.1 was in P + F 2 of C5, and between the same trait and TaCKX11 was in P + F 2 of C3. In seedling roots of various crosses, this trait was mainly negatively correlated with TaCKX5, 8, 10, and NAC2.
Generally, negative correlations between the expression of TaCKX2.1, 2.2.2, and 10 in spikes and TGW, seed number, seed yield, and spike number were correlated with higher yield, and positive correlations were correlated with lower yield. On the other hand, positive correlations between the expression of these genes and root mass determine a higher yield in F 2 . Higher yield in F 2 is also associated with balanced CKX enzyme activity in spikes and seedling roots.
A summary of the regulation of yield-related traits by TaCKX GFMs and NAC2 in the high-yielding progeny of F 2 is presented in Figure 3.

Discussion
Common wheat is a very important cereal crop for feeding the world's population; therefore, continued improvement of the yield of this species is significant. CKX GFMs have already been documented to perform a pivotal role in determining yield-related traits in many plant species, including wheat [4,12]. The genes are tissue-specific; they encode cytokinin oxidase/dehydrogenase, the enzyme that irreversibly degrades cytokinins. We have already characterized the role of TaCKX1 and TaCKX2 in the regulation of yield traits in awnless and owned-spike cultivars [5][6][7]. The range of natural variation in the expression levels of most TaCKX genes among breeding lines and cultivars was very high, indicating the possibility of selecting beneficial genotypes for breeding purposes [8]. Therefore, we were interested in how the expression of these genes is inherited.

The Expression Patterns of Most TaCKX GFMs and TaNAC2-5A Are Mainly Inherited from the Paternal Parent
Comparison of the expression patterns of most of the TaCKX GFMs and yield-related traits between parents and F2 progeny in all reciprocal crosses tested indicated their inheritance from the paternal parent. This rule includes expression patterns in both tissues tested, 7 DAP spikes, and seedling roots, and all TaCKX GFMs and TaNAC2-5A tested were represented in different crosses. The exception was TaCKX5 expressed in 7 DAP spikes, and TaCKX3 expressed in seedling roots, for which the expression level in single crosses was inherited from the maternal parent. Furthermore, high or low yield was predominantly inherited from the paternal parent, and root mass was inherited from both parents or in one reciprocal cross from the maternal parent. We have not found such examples of inheritance in the literature; however, some deviations from parental additivity

Discussion
Common wheat is a very important cereal crop for feeding the world's population; therefore, continued improvement of the yield of this species is significant. CKX GFMs have already been documented to perform a pivotal role in determining yield-related traits in many plant species, including wheat [4,12]. The genes are tissue-specific; they encode cytokinin oxidase/dehydrogenase, the enzyme that irreversibly degrades cytokinins. We have already characterized the role of TaCKX1 and TaCKX2 in the regulation of yield traits in awnless and owned-spike cultivars [5][6][7]. The range of natural variation in the expression levels of most TaCKX genes among breeding lines and cultivars was very high, indicating the possibility of selecting beneficial genotypes for breeding purposes [8]. Therefore, we were interested in how the expression of these genes is inherited.

The Expression Patterns of Most TaCKX GFMs and TaNAC2-5A Are Mainly Inherited from the Paternal Parent
Comparison of the expression patterns of most of the TaCKX GFMs and yield-related traits between parents and F 2 progeny in all reciprocal crosses tested indicated their inheritance from the paternal parent. This rule includes expression patterns in both tissues tested, 7 DAP spikes, and seedling roots, and all TaCKX GFMs and TaNAC2-5A tested were represented in different crosses. The exception was TaCKX5 expressed in 7 DAP spikes, and TaCKX3 expressed in seedling roots, for which the expression level in single crosses was inherited from the maternal parent. Furthermore, high or low yield was predominantly inherited from the paternal parent, and root mass was inherited from both parents or in one reciprocal cross from the maternal parent. We have not found such examples of inheritance in the literature; however, some deviations from parental additivity of expression in polyploid plants were described [2]. An example of such non-additive gene expression takes place when the gene expression level in progeny is higher than that of one parent. The expression level dominance of one parent, also called genomic imprinting, is epigenetic in origin and was investigated primarily at the molecular level in plants and animals [36,37]. The main regulators of gene imprinting are DNA and histone methylation asymmetries between parental genomes. Most of the imprinted genes in the endosperm of grains of different rice cultivars are imprinted across cultivars, and their functions are associated with the regulation of transcription, development, and signaling [39]. Imprinting might affect a single gene or a group of genes. Genes that showed conserved imprinting in cereals have been shown to reveal positive selection and were suggested to regulate seed development in a dose-dependent manner [38]. The only example of a paternally imprinted locus in maize is ded1, which encodes a transcription factor specifically expressed during early embryo development and activates early embryo genes that contribute to grain set and weight [47]. To our knowledge, there are no examples of paternally inherited expression patterns. According to Arabidopsis research, imprinted paternally expressed genes during seed development are mainly related to hypomethylated maternal alleles, repressed by small RNAs or less frequently with transposable elements [48,49]. Contrary to developmental epigenetics, in the case of transgenerational epigenetics, these epigenetic changes do not reset between generations, and this type of inheritance is more related to plants than animals (heritable changes in DNA methylation) [54]. Therefore, we suggest that this paternal inheritance of selected TaCKX GFMs is an effect of transgenerational epigenetic changes, not reset between generations. These heritable epigenetic changes might be effects of DNA methylation, repression of maternal alleles by small RNAs, transposable elements, or, most likely, transcription factors. From our in silico analysis and expression analysis (Iqbal et al., not published yet), several NAC transcription factors appear to strongly regulate the expression of TaCKX GFMs and TaIPT GFMs, influencing yield-related traits [24].

Cooperation of TaCKX GFMs and TaNAC2-5A in the Determination of Yield-Related Traits
The coordinated high or low level of expression of a few groups of genes in the paternal parent positively or negatively regulates higher or lower yield. In two reciprocal crosses, where both TaCKX5 and TACKX9 showed high expression in 7 DAP spikes of the paternal parent, the yield in the F 2 progeny was high and vice versa. In others, the high yield in the F 2 progeny was determined by a low level of expression of TaCKX2.1 and TaCKX11 in spikes of the paternal parent and high levels of their expression in the maternal parent. The level of expression of TaNAC2-5A in the paternal parent and/or F 2 was in opposition to TaCKX5 and 9; however, it was in agreement with TaCKX11 and TaCKX2.1, suggesting their role in the regulation of transcription of these genes. In fact, it was proven by correlation analysis of its expression with yield-related traits [8]. Opposite cooperation of some of the genes in paternal spikes, which resulted in high or low yield in F 2 , has also been observed. The high level of expression of TaCKX2.2.2 and the low level of expression of TaCKX10 predominantly resulted in the low yield in F 2 progeny and vice versa.
Such common rules of gene expression in the paternal parent associated with yield in the paternal parent and F 2 progeny were also observed in the seedling roots. The high level of expression of TaCKX3 and TaCKX8 in the paternal parents of three reciprocal crosses resulted in high yield in the F 2 progeny and vice versa. The expression of TaNAC2-5A in spikes and seedling roots of high-yielding paternal parents and F 2 progeny showed a predominantly low expression level and inversely. These principles of paternal inheritance of selected TaCKX GFMs and TaNAC2-5A expression associated with high yield could be directly involved as molecular markers in high-yielding wheat breeding.

Regulation of Yield-Related Traits by TaCKX GFMs and TaNAC2-5A in the F 2 Generation
The grain number, grain yield, spike number, and TGW were strongly positively correlated with TaCKX2.1 and TaCKX2.2.2 independent of the parent; however, only in the crosses resulted in decreased yield. In contrast, negative correlations were observed between TaCKX2.1, TaCKX2.2.2, and TGW in a reciprocal cross of C7/C8 and TaCKX2.2.2 in a one-way cross (C4), in which the F 2 progeny had a higher yield. All these correlations prove our earlier observations [6,7]. Modified wheat lines with 60% decreased expression of TaCKX2.2.2, and a slight decrease in the TaCKX2.2.1 and 2.1 genes exhibit a significantly higher TGW and slightly increased yield [7]. Interestingly, this result was observed in cultivars and breeding lines that represent awnless spikes. In the owned-spike cultivar, silencing of the TaCKX2 genes co-expressed with other TaCKX resulted in decreased yield; however, TGW was at the same level as in non-silent plants [5]. Furthermore, a strong feedback mechanism for regulation of the expression of TaCKX2 and TaCKX1 genes was observed in both awnless and owned-spike cultivars [5][6][7]. Silencing of TaCKX2 genes upregulated the expression of TaCKX1 and vice versa. This feedback mechanism could explain the observed positive correlations of the TaCKX2 genes with yield-related traits in low-yielding F 2 progeny and negative correlations in high-yielding F 2 progeny. A similar mechanism is visible when we analyze individual traits in high-yielding F 2 , such as grain number, grain yield, and spike number. These traits are promoted by up-regulated in 7 DAP spikes TaCKX2.1 and down-regulated TaCKX1. Silencing of HvCKX1 in barley, which is an ortholog of TaCKX1, decreased CKX enzyme activity and led to increased seedling root mass and higher plant productivity [55]; however, knock-out of this gene caused a significant decrease in CKX enzyme activity but no changes in grain yield were observed [56]. These differences might be explained by differences in the level of decreased gene expression, which variously coordinate the expression of other genes, regulate phytohormone levels, and determine particular phenotypes, as was already documented in wheat [5,7].
Grain number was also negatively correlated with TaCKX1 and TaCKX5, and grain yield was negatively correlated with TaCKX1, predominantly for M and F 2 , which were characterized by decreased grain number and lower yield. This observation is also in agreement with previous research. The silencing of TaCKX1 caused an increase in spike number and grain number but a decrease in TGW because this trait is opposite to grain number [6]. The low-yield progeny of F 2 showed positive correlations between the expression of TaCKX11, 3, 5, and 8, and NAC2 in the seedling roots and the grain number, the spike number and the grain yield; however, the higher-yield progeny of F 2 displayed a negative correlation between TaCKX1 and these yield-related traits in the seedling roots of some crosses. In summary, high-yielding F 2 was the result of upregulation of TaCKX2.1 in spikes and downregulation of TaCKX1 in seedling roots. As documented earlier, TaCKX11, 5, 8 and TaNAC2-5A are expressed in all organs, and their expression is correlated with the expression of spike-specific TaCKX2 and TaCKX1 [8,58].
Rice OsCKX11 is an orthologue of wheat TaCKX11 and is highly expressed in the roots, leaves, and panicles. The gene was shown to coordinate the simultaneous regulation of leaf senescence and grain number by the relationship of source and sink [60]. Since TaCKX11 is expressed in seedling roots and highly expressed in leaves, inflorescences, and 0, 7, and 14 DAP spikes, it could perform a similar function. This is partly proven by silencing of the TaCKX2 genes in awnless spikes of cv. Kontesa, which resulted in significant upregulation of TaCKX11 and growth of TGW, and chlorophyll content in flag leaves [7]. In contrast, TaCKX11 is significantly negatively regulated by TaCKX1, resulting in a higher spike number and grain number [6]. Its orthologue in rice, OsCKX11, was found to regulate leaf senescence and grain number by the coordinated source and sink relationship [60].
Based on a summary of the regulation of yield-related traits in high-yielding F 2 , it is possible to identify singular genes or groups of genes that are up-or down-regulated in 7 DAP spike or seedling roots and specifically regulate yield-related traits. Upregulated in spikes TaCKX2.1 and downregulated in seedling roots TaCKX1 were found to determine grain number, grain yield, spike number, and spike length. Furthermore, upregulated in 7 DAP spikes TaCKX10 and downregulated TaNAC2-5A, together with others, depending on cross, control spike length, semi-empty spikes, root mass, and increased grain yield. As discussed above, high TGW is in contrary to high grain number and partly grain yield and was strongly determined by downregulated TaCKX2.2.2 together with TaCKX2.1 in 7 DAP spikes and upregulated TaCKX10 and NAC2 in seedling roots. The upregulated in seedling roots TaNAC2-5A participates in the determination of TGW and plant height, and the downregulation of TaNAC2-5A in seedling roots controls the development of semi-empty spikes and root mass.
In previous research, TaNAC2-5A has been documented as a gene encoding a nitrateinducible wheat transcription factor. Overexpression of the gene improved root growth, grain yield, and grain nitrate concentration [25]. This is in agreement with our observations of growth of TGW but not enhanced roots. The increase was argued to be the consequence of regulation of nitrate concentration and its remobilization in developing grains by direct binding of the TaNAC2-5A protein to the promoter of the nitrate transporter, TaNRT2.5-3A and positive regulation of its expression [61]. The expression of TaNAC2-5A is coregulated by expressed in 7 DAP spikes TaCKX2 genes and expressed in 7 DAP spikes and seedling roots TaCKX1 gene [5,7]. Independent of awnless or awned-spike genotype, downregulation of TaCKX2 genes by RNAi significantly increased TaNAC2-5A expression, resulting in higher chlorophyll content in flag leaves and delayed leaf senescence. As discussed above, the strong feedback mechanism between the TaCKX2 and TaCKX1 genes implies that downregulation of TaCKX1 resulted in opposite results. Similar to our observations in wheat, an ortholog of TaNAC2-5A in rice, OsNAC2, was described as a negative regulator of crown root number and root length [62]. Its expression was positively correlated with cytokinin synthesis genes, OsIPT3, 5, the gene determining the formation of active cytokinins, OsLOG3, and negatively correlated with OsCKX4 and 5. The authors concluded that OsNAC2 stimulated cytokinin accumulation by suppressing CKX expression and stimulating IPT expression by binding the OsNAC protein to the promoters of these genes. Therefore, OsNAC2 functions as an integrator of cytokinin and auxin signals that regulate root growth. In our experiments, orthologous to OsCKX4, TaCKX4 was not tested due to its weak expression in roots. However, the up-regulated expression of highly specific in seedling roots TaCKX3 and TaCKX8 [8,58] was antagonistically regulated by TaNAC2-5A in these organs, positively influencing seedling growth. Furthermore, our in silico analysis of TaNACs with TaIPTs and TaCKXs showed that the same NAC proteins might join promotor sites of cytokinin synthesis and cytokinin degradation genes in wheat (Iqbal et al., not published yet).
Ten seeds of each genotype germinated in Petri dishes for five days at room temperature in the dark. Six out of ten seedlings from each Petri dish were replanted in pots with soil. The plants were grown in a growth chamber under controlled environmental conditions with 20 • C day/18 • C night temperatures and a 16 h light/8 h dark photoperiod. The light intensity was 350 µmol s −1 ·m −2 . The plants were irrigated three times a week and fertilized once a week with Florovit according to the manufacturer's instructions.
The following tissue samples in three biological replicates were collected: 5-day-old seedling roots, which were cut 0.5 cm from the root base before replanting in the pots, and first 7 DAP) spikes from the same plants grown in the growth chamber. All of these samples were collected at 9:00 a.m. The collected material was frozen in liquid nitrogen and kept at −80 • C until use.

Cross-Breeding
The maternal plant was deprived of its own anthers so that it would not self-fertilize, then pollinated by transferring three anthers from the paternal plant for each ovary of the maternal parent plant and placed in an isolator. The seeds were harvested.

RNA Extraction and cDNA Synthesis
Total RNA from 7 DAP spikes and roots from 5-day-old seedlings was extracted using TRI Reagent (Invitrogen, Lithuana) according to the manufacturer's protocol. The concentration and purity of the isolated RNA were determined using a NanoDrop spectrophotometer (NanoDrop ND-1000, Thermi Fisher Scientific, Wilmington, DE, USA), and the integrity was checked on 1.5% (w/v) agarose gels. To remove residual DNA, RNA samples were treated with DNase I (Thermo Fisher Scientific, Lithuana). Each time, 1 µg of good quality RNA was used for cDNA synthesis using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Lithuana) following the manufacturer's instructions. The cDNA was diluted 20 times prior to use in the RT-qPCR assays.

Quantitative RT-qPCR
RT-qPCR assays were performed for 10 genes: TaCKX1, TaCKX2.1, TaCKX2.2.2, TaCKX3, TaCKX5, TaCKX8, TaCKX9, TaCKX10, TaCKX11, and TaNAC2-5A. The sequences of the primers for each gene are shown in Table S2. All real-time reactions were performed on a Rotor-Gene Q (QIAGEN Hilden, Germany) thermal cycler using 1× HOT FIREPol EvaGreen qPCR Mix Plus (Solis BioDyne, Estonia), 0.2 µM of each primer and 4 µL of cDNA in a total volume of 10 µL. Each reaction was carried out in three biological and three technical replicates in the following temperature profile: initial denaturation and polymerase activation of 95 • C-12 min (95 • C-25 s, 62 • C-25 s, 72 • C-25 s) × 45 cycles, 72 • C-5 min, with melting curve at 72-99 • C 5 s per step. The expression of TaCKX genes was calculated according to the two standard curve method using ADP-ribosylation factor (Ref 2) as a normalizer. The relative expression for each TaCKX GFM and TaNAC2-5A was calculated in relation to the control female parents, set as 1.00.

Analysis of CKX Activity
CKX enzyme activity was performed in the same samples subjected to TaCKX gene expression analysis according to the procedure developed by Frebort et al. [63] and optimized for wheat tissues. The plant material was powdered with liquid nitrogen using a hand mortar and extracted with a 3-fold excess (v/w) of 0.2 M Tris-HCl buffer, pH 8.0, containing 1 mM phenylmethylsulfonyl fluoride (PMSF) and 0.3% Triton X-100 ((St. Louis, MO, USA). Plant samples were incubated in a reaction mixture consisting of 100 mM McIlvaine buffer, 0.25 mM of the electron acceptor dichlorophenolindophenol and 0.1 mM of substrate (N6isopentenyl adenine). The volume of the enzyme sample used for the assay was adjusted based on the enzyme activity. The incubation temperature was 37 • C for 1-16 h. After incubation, the reaction was stopped by adding 0.3 mL of 40% trichloroacetic acid (TCA) and 0.2 mL of 2% 4-aminophenol (PAF). The product concentration was determined by scanning the absorption spectrum from 230 nm to 550 nm. The total protein concentration was estimated based on the standard curve of bovine serum albumin (BSA) according to the Bradford procedure [64].

Measurement of Yield-Related Traits
Morphometric measurement of yield-related traits of selected genotypes was performed. The described traits were plant height, spike number, semi-empty spike number, tiller number, spike length, grain yield, grain number, TGW, and 5-day seedling root weight.

Statistical Analysis
Statistical analysis was performed using Statistica 13 software (StatSoft). The normality of the data distribution was tested using the Shapiro-Wilk test. The significance of the changes was analyzed using ANOVA variance analysis and post hoc tests. The correlation coefficients were determined using parametric correlation matrices (Pearson's test) or a nonparametric correlation (Spearman's test).

Conclusions
We indicate, for the first time, that the pattern of expression of selected TaCKX GFMs and TaNAC2-5A, and grain yield in wheat, is paternally inherited by the F 2 generation. Pater-origin transmission of gene expression levels sheds new light on the method of parent selection and crossing to obtain high-yielding phenotypes. We also showed which genes cooperate together by upregulation or downregulation and which function in the opposite manner in establishing yield-related traits. This knowledge can be applied to select the desirable phenotype in F 2 . For example, a high-yielding paternal parent with downregulated, compared to the maternal parent, expression of TaCKX2.1 and TaCKX11 in 7 DAP spikes and upregulated expression of TaCKX3 and TaCKX8 and downregulated TaNAC2-5A in seedling roots is expected to transmit this pattern of expression to F 2 , which will result in a high yield. The main problem is the antagonistic expression patterns of genes for some important yield-related traits, such as grain number, grain yield, and spike number, to TGW, which is the result of the feedback mechanism of the regulation of expression between TaCKX1 and TaCKX2 genes and others. The expression analysis of TaNAC2-5A and the in silico analysis of TaNAC GFMs revealed that the encoded proteins participate in the regulation of transcription of selected TaCKX genes responsible for cytokinin degradation and TaIPT genes responsible for cytokinin biosynthesis. Therefore, TaNACs are important additional regulators of yield-related traits in wheat, which should be taken into consideration in wheat breeding.